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Abstract. Wc discuss the constraints imposed on the nonhnear evolution of the 
Large Scale Structure (LSS) of the universe by galilean invariance, the symmetry 
relevant on subhorizon scales. Using Ward identities associated to the invariance, 
wc derive fully nonlinear consistency relations between statistical correlators of the 
density and velocity perturbations, such as the power spectrum and the bispectrum. 
These relations are valid up to 0{f%]^) corrections. We then show that most of the 
semi-analytic methods proposed so far to resum the perturbative expansion of the 
LSS dynamics fail to fulfill the constraints imposed by galilean invariance, and are 
therefore susceptible to non-physical infrared effects. Finally, we identify and discuss 
a nonperturbative semi-analytical scheme which is manifestly galilean invariant at any 
order of its expansion. 
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1. Introduction 

The Large Scale Structure of the Universe (LSS) provides an unique arena to test our 
understanding of the evolution of the Universe in recent cosmological epochs. The main 
obstacle towards an accurate connection between theory and observations is represented 
by the many different sources of nonlinearity: the nonlinear evolution of the dark matter 
(DM) fluid, redshift space distortions, halo and galaxy bias, etc. 

Cosmological perturbation theory has received considerable attention in the recent 
literature as a promising tool to approach the nonlinear scales [|] Being based on analytic, 
or semi-analytic, techniques, its computational times are in general greatly reduced with 
respect to those for N-body simulations. Moreover, although the basic formalism is 
derived for an Einstein-de Sitter cosmology, its extension to ACDM is straightforward, 
and also its formulations in non-standard cosmologies are feasible and under control. 

The main limitation of standard perturbation theory (PT) [Ij is in the range of 
scales in which it can provide accurate results. As it was shown in [2], percent accuracy 
in the DM power spectrum (PS) can be attained for wave vectors k ~ 0.2hMpc~^ for 
redshift z > 1, rapidly degrading at high k and smaller z. Therefore, it cannot provide 
accurate predictions for the PS in the baryon acoustic oscillation (BAO) range of scales 
(0.05hMpc~^ ~ ~ 0.25hMpc~^) at low redshifts. However, starting from the works 
of Crocce and Scoccimarro [31 H], it has been realized that the PT expansion can be 
reorganized very efficiently, so that all-order PT contributions can be taken into account 
at lowest orders in the new expansion schemes [SI El El El El [IHl Ell [121 [121 [IH [U] • 

Progress along these lines of research has been quite impressive: independent 
approaches are now able to produce, in a 0(1 min.) time, DM PS which agree with 
those from high precision N-body simulations (and among themselves) at the percent 
level in the full BAO range of scales down to 2; = [TH |T71 [HI [19], and to a few percent 
at smaller scales, up to A; ~ 0.8hMpc~^ |17j . 

The purpose of this paper is to fully exploit the symmetry of the system, namely 
galilean invariance (GI), as a portal towards the nonlinear scales. The use of symmetry 
principles as tools to constrain the perturbative and nonperturbative sectors is common 
practice in quantum field theory. The consequences of the fundamental symmetries are 
enforced by Ward identities and consistency conditions, and provide powerful constraints 
on the structure of the fully renormalized (or effective) field theory. These constraints 
must be satisfied, besides perturbation theory, by any other viable approximation 
scheme. 

GI in the context of the LSS was discussed in |2(Jj. Since a galillean transformation 
(GT) corresponds to a velocity perturbation of infinite wavelength, the link between GI 
and the infrared (IR) sector of the PT expansion was elucidated. In particular, the fact 
that leading IR divergences (emerging for a scale-free PS P ~ A;" with n < — 1) cancel 
out at any order in PT [21j, was clearly understood as a consequence of GI. This fact 

I In this paper, we will consider only the first source of nonlinearity mentioned above, namely the 
treatment of nonlinearites in the evolution of the pure DM fluid in real space. 
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was then used to identify which modifications of the full Euler-Poisson system respect 
GI and are therefore free of IR divergences, selecting the Zel'dovich dynamics as the 
only viable approximation. 

In this paper we build upon this previous work, with a different focus. Namely, 
we stick to the exact dynamics, described in terms of the Vlasov equation and the 
Euler-Poisson system, and: 

• we derive the Ward identities, or consistency relations, enforced by GI on the 
correlators, at the fully nonlinear level. An example of it is the relation between 
the fully non-linear bispectrum and the fully non-linear PS of the density contrast 
S given by eq. ( 74 ) : 



/c-s>0 



lim Bs{k,q, |q + k|;r, t',t") 

q-kD+{T')-D+{r") Qf^^Q^^pH^, D+ {tq) 

P D+{t) k'^T{k) D+(t) 

0{k',fL) , 



Ps{k;T,T)Ps{q;T',T") 



(1) 

(see Sects. [2] and |5] for the definitions of the various quantities). While it is 
straightforward to verify this relation at the tree level. Ward identities guarantee 
that this relation holds at the full non-linear level. This is completely analogous to 
the non-renormalization of the electric charge in QED through the nonperturbative 
Zi = Z2 identity. Finally, we stress that this relation also accounts for the present 
dark energy; 

• we show that the nonperturbative semi analytic approaches mentioned above are 
generally not well-behaved from the GI point of view, and that, indeed, most of 
them are plagued by spurious IR effects; 

• we identify and discuss a galilean invariant PT resummation scheme, the 'eikonal 
renormalized perturbation theory' (eRPT), introduced in [T7]. 

The paper is organized as follows. In Sect. [2} we introduce GI in physical and 
comoving coordinates, where it assumes the form of a time-dependent velocity boost. 
Then we show that both the Vlasov equation and the continuity and Euler equations 
derived from it are GI. In Sect, [s] we review the path- integral approach of |S] as an 
useful formalism to discuss symmetry properties at the fully nonperturbative level. In 
particular, we show how to extend this formalism from the center of mass (c.o.m.) frame 
to any frame related to it by a GT. Then, in Sect. |4]we show how to derive the Ward 
identities of GI, and we focus on those linking the two- and three-point correlators. In 
Sect. [5| we discuss how this relation is affected by an initial nongaussianity, leading to 
the result written above. In Sect. [6] we discuss the link between GI and the IR sector 
of the loop corrections. We show how a GT can be implemented in a diagrammatic 
language and, finally, introduce a procedure to check if a given approximation scheme 
respects the GI constraints on the IR sectors and is therefore free from spurious IR 
effects. The link between GI and IR sensitivity is further elaborated in Sect. [7j In 
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Sect. IHlwe consider the PT resummation schemes and show that most of them are not 
free from spurious IR effects. On the other hand we show that the eRPT [17] respects 
both our test procedure and the Ward identities at any order in the expansion. Finally, 
in Sect. [9| we discuss the implications of our results. 



2. Galilean invariance of the tree level action 



A GT in physical coordinates reads: 
R' = R - w t 

V' = V - w . (2) 

Comoving coordinates and conformal time are defined, respectively, as 
R 

X = — , 
a 

dT=-, (3) 
ct 

where a(r) is the scale factor. Moreover, the peculiar velocity and its canonically 
conjugated momentum are v = dx/dr and p = am\. The transformations (|2| then 
take the form 

x' = X — w T , 

v' = V — w T , (4) 

where 

T(r) = ^ TdrVr'), 
a[V Jo 

■ , , dT f 1 da\ 

r(r) = - = i-«r («^-^). (6) 

that is, in comoving coordinates and conformal time, GT can be seen as a boost by 
a time- dependent velocity w(r) = wT(r). This makes the consequences of GT on 
multi-point correlators non-trivial. 

We note that, under the transformation Q, the peculiar momentum transforms as 
p' = p — arriTw. 

The particle distribution function is a GT-invariant scalar: /'(x', p', r) = 
/(x, p, r) = /(x' + w T, p' + arriTw, r), and the total time-derivative operator, 

d d ■ d ■ d 

is GT-invariant. Therefore, the Vlasov equation, namely, 

^/(x,P,r) = ^/'(x',p',r), (7) 



is GT-invariant. Indeed 
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(I- + K - w'T)^ + - ^{amf)w')^ 
^dr ^ ' dx' ^ dr^ ' ' dp^ 

. . d d ■ ■ d 

+ ^"^0^ + -^i(^^T)w')-^)f{^' + wT, p' + amTw,T) 



(8) 



The Vlasov equation can be cast in a series of equations for the moments of the 
distribution function [T] 



J rf^p/ (x, p,t)=p (x, r) = po (t) [1 + 5 (x, r)] 

/ rf3p^/(x,p,r) 
J am 



p [X, T) Vi [X, T) 



(9) 



The series can be consistently truncated to the continuity and Euler equations by setting 
aij = (this truncation being known as the "single stream approximation"): 

' d 



Or 
d_ 



+ V ■ V 
+ V ■ V 



6 



[1 + 5) Vv 



V = - V0 , 



(10) 



where (f) is the gravitational potential, which, on super-horizon scales, obeys the Poisson 
equation 

_o , 3, 



-1~L ^rn5 



11 



where f2m(fl) is the time dependent matter density parameter. Throughout this paper 
we will assume a ACDM background cosmology, or, more generally, a dynamical dark 
energy scenario in which the dark energy fluctuations can be neglected on all scales of 
interest. 



The expressions (10) have been written in such a way that both LHS's and RHS's 
are invariant under a GT. We note that on the LHS's both the linear and nonlinear 
terms transform nontrivially, but that the two nontrivial pieces in the transformation 
cancel each other, so that the total operator is invariant (the cancellation is inherited 
from that of the Vlasov equation shown in ([s])). We also note that the gravitational 
potential transforms nontrivially under a GT |20j. 

This system can be reformulated with a compact notation [31 HJ [5] by Fourier 
transforming 



5(x,r) 



d'p e'"-'' 5 (k, r) , 



and analogously for v, by assuming no vorticity, 

k 



v(k,r) = -z-^^(k,r) 



(12) 



(13) 
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with fc^ = k ■ k, and by introducing the doublet 



V2 (k,r7) 



5(k,r/) 



7] = In 



D+{Tin) 



(14) 



In the definitions above, we have introduced the new time variable, t], in terms 
of the linear growth factor and we have also used the linear growth function 
f{r]) = l/n dr]/dT. 

With this notation, eqs. ( 10 ) can be reformulated as 

{Sabdr, + ^ab) Vb {K V) = / ^^^P d'^'ilabc (k, -p, -q) (p, T]) (q, T]) , (15) 



where 



17 



1 -1 

"2 p 2 p 



(16) 



and where the only non vanishing components of the vertex function are 

7121 (k, p, q) = (k + p + q) a (p, q) 7i2i (k, p, q) = 7112 (k, q, p) 
7222 (k, p, q) = 5d (k + p + q) /3 (p, q) , 
with 5d being the Dirac 5— function, and 

(p + q) ■ p 



(17) 



a (p, q) 



/3(p,q) 



ip + q) p q 

2p2g2 



Under a GT the field transforms as 

(k, r/) e^''--^(''Va (k, r/) + z k ■ w e'^^.T (r^) 5d (k) 5,2 ■ 



(19) 



The second piece must be retained in the GT of the right hand side of ( 15 ), as the vertex 
functions have a pole when their second or third argument vanishes. Their presence is 
actually crucial to ensure the covariance of (15). Indeed, under a GT, eq. (15) changes 
into 

{Sabdr^ + ^ab) [c^*^'^ ^("Vfe (k, r/)] = j d3pd3q^^^^(k, -p, -q)e^(P+'i)-^('') x 

[Pb (p, ??) + «P ■ we'''^dr,T6D (p) ^62] [fc (q, r?) + iq ■ we'^^d^TdD (q) 5c2] • 

(20) 

There are four terms at right hand side. The one obtained from the two second pieces 
in the square parenthesis vanishes, since the various 5— functions force k = 0, and the 
vertex function vanishes when its first argument vanishes, due to translation invariance. 
The two "mix-terms" instead add up to cancel the derivative of the phase on the left 
hand side. Indeed, using 

1 p ■ qp ■ w 



P-w6d (p) labc (k, p, q) 



-6b2SacSD (k + p + q) (5b (p) , (21) 
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eq. (20) rewrites 



^ik-w T(r?) 



e^^ d pd q'jabc (k, -p, -q) (p, r^) (q, 



+ a^T(^4k,r/) / d^p 



ip ■ w p ■ k 

p2 



(22) 

We indeed see that the last term at the LHS cancel against the last term at the RHS, 



and eq. (15) is recovered also in the new frame. We stress that the two terms that 



cancel each other come from the transformation of the linear and the nonlinear part of 



(15). 



As done in [5] we introduce the action 

dridri' d^kxa (-k, ri) g-^{r], r]') (pb (k, r/) 



^ dr]d^kd^p d^qe'' '^abc (-k, -p, -q) Xa (k, t]) (fb (p, v) fc (q, v) ,(23) 
where 

9abM) = Sd{V - V') iSabd,' + Qab) , (24) 

is the inverse linear propagator, which, inverted imposing casual initial conditions, gives 
the linear propagator, 

gabiv)={Bab + e-'^''^Aab) 9(7]), 

with 9{x) the Heaviside step function, and 

5 I 3 2 / ' 5 I 3 



B 



(25) 



(26) 



In order to obtain the above explicit expression, we have approximated fim//^ 



in (16), see [U HO] for a discussion of this approximation and for an assessment of 
its numerical validity at the non-linear level. This approximation allows for a simple 
expression of the linear propagator, but it is not used in the reminder of the paper, and 
therefore it by no means affects the following discussion. 



Extremizing the action in eq. (|23|) with respect to Xa, gives the equation of motion 

(27) 



(15). To enforce invariance under GT, we impose that, under a GT 
k,r/)->e^'^--^('')Xa(k,r^) , 



Xa 



SO that, under (19) and (27) 



S ^ S 



J d7]d^kxa (-k, v) {Sa2d^n + ^a2) ^k ■ we-^S.T^z) (k) 



(2J 



which is invariant up to a term that gives a vanishing contribution to the equation of 



motion (15). The invariance is the result of the same cancellation between the linear 



and the nonlinear term that we have already observed in (22), and that we have also 



commented after eq. (11). 
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3. Path integral formulation and frame fixing 



Ref. [5] provided a path integral formulation of the system (15). The formulation is 
characterized by the generating functional 

Z [Ja, K,- Vo] = j V^aVxb exp I - ^ y rf^k xa (-k, 0) {k) Xb (k, 0) + iS 

+ i j dr^rf^k [J„ (-k, ri) (k, 7]) + Ka (-k, r/) Xa (k, r])] + iS^ [Vq] 



(29) 



where S is the action (23). Correlators with n ifa fields and m Xb ones are obtained, 
as usual, by taking n derivatives of the path integral with respect to the source iJa and 
m derivatives with respect to iKfy. In particular, second derivatives give the nonlinear 
two-point correlator, i.e. the PS, 
5W 

i{ipa (k, r^) (k', V))vo = t6D{M)Pab{k- 7], v') , (30) 



5Jai-'k,r])6Jbi-k',T]') j^^^^^q 
and the nonlinear propagator 
6W 



t{ipa (k, r/) Xb (k', v'))vo = -SD{\^+\<.')Gab{k; v, v') , (31) 



Ja=Ka=0 



where W = —i log Z. 

The expression in eq. (29) holds for gaussian initial conditions (at the "time" ?7 = 0, 
chosen when all the relevant modes are linear), encoded by the linear PS 

P'^,ik) ^ P,V^; 0, 0) = uau, P° (k) , (32) 



where in the last equality Ua 



selects the growing mode. Nongaussian initial 



conditions can be taken into account by adding trinliear, quadrilinear, and so on, terms 
in Xa for the initial bispectrum, trispectrum, etc., respectively, see eq. (62). 

In (29) we introduced also a term linear in Xa, which was absent in [5], the "frame- 
fixing" term 



5ff [Vol 



d^kdr]dr]'xa (-k, r]) gj {r], t]') (ph (k, r]', Vq) 



(33) 



where 



ipa (k, T], Vo) = -I (5,26-^ k-YoSo (k) d,T. (34) 

Adding the term Ss to the action selects the inertial frame in which the average physical 
velocity of particles is given by Vq, and, correspondingly, the average peculiar velocity 
is VoT, see eq. nAn. Indeed, we first notice that 



(v(x)) 



Vol. 



d xv(x) 



^^Hfe^ Jd'k6n{k) ^^^(k) 



(35) 



where we have used eqs. (13) and (14) and "Vol." stands for the total volume on which 
the average is taken. Setting (p2(k) = (p2 (k, r], Vq), gives (v(x)) = VqT (where we have 
used the relation T-LfdrfT = T). 



Galilean invariance and the consistency relation for the nonlinear squeezed bispectrum of large scale structur 



Then, we show that the inclusion of the hnear term in Xa, eq. (33), at the 



exponent of the path integral forces the field V5a(k, rj) to have a non- vanishing expectation 
value, given by eq. (34). Indeed, when Sg = 0, translational invariance of the action 



and of the initial conditions ensure that the expectation values of the field vanish. 
Translational invariance of the dynamics is encoded in the property of the linear vertex 

limfc^o 7a6c(k, q, p) = 0. 

To see how a non-vanishing expectation value emerges when Sg ^ we transform 



it away by making a change of integration variables in (29) corresponding to eqs. (19) 
and (27) with w = — Vq. It gives 

Z [J„, Ki,; Vo] = e^/'^''^'''^"(-''''')^'"('^''''V°) Z [j, e'^^''^^''), e'^^°^^'^^;Vo = O] , (36) 

where the generating functional on the RHS is obtained from eq. ( 29 ) by setting 5*^ = 



and by multiplying the sources by the GT phases. The field expectation value is obtained 
by taking the first derivative with respect to the source, and then setting the sources to 
zero, 

dlogZ[Ja,K,;Vo] 



(37) 



Ja=Ka=0 



dJa{-'k,ri) 

The only non- vanishing contribution to the derivative of the RHS of (36) comes from 
the factorized exponential, which gives precisely 



{^a{Kv))vo = Va (k,r/,Vo) 



(3J 



By taking the appropriate derivatives of eq. (36) with respect to the sources, we can 



obtain the transformation law for any correlator of (p and x fields between the center of 
mass frame (Vo = 0) and a generic one with Vq 7^ 0. It is given by 

{fa^ (qi, ^1) • • • ^ar, {^.n, Vn) Xa„+i i<ln+U Vn+l) ■ ■ ■ Xa„+™ {^.n+m, Vn+m))vo 
X ifai (qi, ^1) • • • fa„ (qn, Vn) Xa„+i (qn+1 

(39) 

where we have taken all momenta 7^ 0. Since translational invariance requires the 
correlator to be proportional to Sd (^^qi), we see that equal-time correlators, such as 



the PS from eq. (30) with rj = rj', are invariant under a GT. 



4. Ward identities 



In Sect. 2 we have seen that the the covariance of the equations of motions, eq. (15) 



under GT, emerges as a consequence of a relation between the tree-level (linear) inverse 
propagator and the tree-level vertex in the limit in which one of its y^a-ends corresponds 



to a velocity field in the long wavelength limit (see eqs. (20), (21)). In this section, we 



investigate how the galilean invariance of the dynamics (and of the initial conditions) 
manifests itself at the fully nonlinear level. We will find relations between correlators 
holding nonpertubatively, and therefore valid beyond the range of scales in which usual 
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perturbative (standard and renormalized) methods can be apphed. Using a field- 
theoretical language, these "consistency relations" are the Ward identities of GT. 

Let us consider a GT of infinitesimal parameter w. Under this transformation, the 



fields (fa and Xa transform as the infinitesimal versions of (19) and (27), respectively 



5Lpa (k, 1]) =ik w [T [t]) ifa (k, r/) + e ''S^T (r/) 5d (k) ^^a] 
^Xa (k, ri) =ik-wT {r]) Xa (k, rf) , 
while the generating functional Z remains invariant: 



(40) 



5Z = j V(fiaT>Xbexp{...} j drid^kik 



W X 



Ja (-k, 7]) (k, ri) + Ka ("k, 7]) Xa (k. If]) 



+ 



J2 (-k, 7]) + Xa (-k, 7]) {Sa2dr, + ^^02) 



e-'^d,T{7])6D{k) \=0, (41) 



where the ellipses in the exponent denotes the terms present in Z = J T>ipa T^Xh exp {•••}■ 
The IPI Green functions are obtained from functional derivatives of the effective 
action F, related ioW = —i log Z by a Legendre transform: 



r [^a, xA = W [Ja, Kb] - J d7]d^k [Ja (-k, T]) ^a {V, k) + Ka (-k. 



where the 'classical' fields in presence of sources are 

6W 



and 



'fa (KV) 



Jaik,ri) 



SJn 



-k,7]) 
5T 



, XaiKv) 



6W 



5Ka 



Kaik,ri) 



-k,7]) ' 
5T 



r/)Xa(r/, k)] , (42) 



(43) 



6!fa{-k,7]) ' ""^"'''^ 6Xa i-k,7])' ^^^^ 

In the following, we will omit the overtilde on the classical fields but it should be clear 
that they are (source-dependent) expectation values, not be confused with the stochastic 



fields we have considered up to now. Starting from (41), one can take out of the path 



integral all the field independent quantities, and express the sources through (44). One 
IS left with[§] 



j dr]d^k (k ■ w) |r {rj) 



5T 



+ 



6T 



Xa 



5fa {KV) 
k, 7]) {5a2 + ^a2 



Va (k,?]) + 



5T 



^Xa (k,?7) 
e-'^dr,T{7]) 6d (k) 



Xa (k, rj) 
= 0. 



(45) 



_6ip2 {k,7]) 

Different one-particle irreducible (IPI) Green functions can be related to each other 
by taking functional derivatives of this expression, and then setting the fields to zero. An 



analogous expression, in terms of W and the sources, can be obtained from (41) using 



We note that we have disregarded the infinitesimal transformation of in cq. (41). Using eq. (36 1 



we can see that the effect of this term is that, in a generic frame, all the J and K appearing in equations 
(41 1 to (44 1 should be multiplied by the phase e^^^°'^ . The expression (45) and the following ones are 



then unchanged. 
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(43), which will generate relations between connected correlators. The physical contents 
of the latter are of course identical to that of the equations for IPI functions. We discuss 
the relations obtained in terms of IPI functions, since they show more directly how the 
tree-level relation between different terms in the action (the inverse propagator and the 
vertex) are generalized at the nonperturbative level. 

As a first example, we take the S'^/S(j)bSxc derivative of eq. (45), and integrate by 
parts in rj the term coming from the second line. We get 

5^T 



[p ■ w 5b (r/ - V) + q ■ w 5d (?7 - r]")) 



Va=Xb=0 



- / d^kk ■ w 6d (k) e"" 



6^T 



6(pb (p, V') ^Xc (q, V") (k, T]) 
Multiplying by and integrating in t], we get 



. (46) 



'Pa=Xb = 



+ J dr]d'kk-w5Dik)TZ''i-p-Kp,^;v'\v',v) = 0, (47) 
where we have used the definition of the inverse (fully nonlinear) propagator [S] 

= Sniq + P) {gba'iv", v') - ^Up; v", v')) , 

5D(q + p)G,-;(p;V',V), (48) 



(p, v') ^Xb (q, v'' 



and of the full trilinear vertex 
5^T 



SXa (k, 1]) 6ipb (p, r]') 6(pc (q, r]") 



= Soik + q + p)T^'^f{-q - p, p, q; r/, r/', r]") 



'Pa=Xb=0 



-2Sd{v - v')Sd{v - V") labdK p, q) + 0(1 - loop) . (49) 



Eq. ( 47 ) is one example of the analogous of Ward identities in quantum field theory. It 
enforces the constraint on fully renormalized {i.e. nonlinear) quantities (in this case, the 
inverse propagator, and the trilinear Xff vertex) coming from the underlying symmetry 
of the theory, in this case galilean invariance. Using the explicit expression for the linear 
propagator in eq. (24) and eq. (21) for the linear vertex, one can immediately check the 
identity (47) at tree-level. 

At one loop, the relation (47) acquires the diagrammatic form shown in Figure [l| 
where continuous lines with an open square represent linear PS's and continuous-dotted 
lines represent propagators, see [5] for details. In the A; — t- limit, and using again the 
property (21 ), we get that the contributions from the upper tree-level vertex in the three 

k- 1 



diagrams at the RHS sum up to give 



A;2 



A;2 



kp 

A;2 



(50) 



where 1 is the loop momentum. Since the sum is 1- independent, it factorizes from the 
loop integral, which gives exactly the 1-loop expression for T,cb{p', v"^ v') ^^^d, multiplying 
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eq. (50| by k ■ w SdO^), and integrating in d^k reproduces the p ■ w (e^" — e''') factor at 
the LHS of eq. gg. 

Taking the S'^/Sxb^Xc double derivative of (45), and performing similar 
manipulations as those used in obtaining eq. (47), we get another relation between 
two- and three-point IPI functions. 



+ j dr7rf'kk-w5,,(k)r^,Y(-k-q,q,k;r/",V,r/) = 0, (51) 



where [5] 



5^1 



^Xa (q,^) ^Xb (p,V) 



i5d{(\+V) {P\q)uaUb6Div)5Div') + ^abiq;V,v')) ,(52) 



'Pa=Xb=0 



and 



6^ 

SXa (k, T]) 6xb (q, v') (p, v") 



= SdO^ + q + p)Tabci-^ - p, q, p; V, v', v'' 



Va=Xb=0 



(53) 



Notice that the trilinear vertex T-^^f is absent at tree-level, but, as 



6^r 



fa=Xb = 



Sxa (k, T]) 6xb (q, v') ^Xc (p, n 

it is generated at 1-loop and higher orders 



= 5D(k + q + p)r^6^c^(-q - p, q, p; r], r]' , r]" 



(54) 



Ward identities (47) and (51) can be combined to derive a relation between the 
nonlinear bispectrum and the nonlinear power spectrum. Indeed, the fully nonlinear 
bispectrum 

{ipa{k;ri)ipb{q,]ri')ipc{p;ri")) = 60^^ + q + p)Babc{k,q,p;r],T]' ,r]") , (55) 

can be exactly expressed in terms of the full trilinear IPI vertices, the full PS and the 
full propagator, as follows 



Babc{k,q,p;r],ri',r]" 



ds ds ds" 



- {Ga,{k- r/, s)PbM\ V', s')Pcf{p; v", ^")rr/(k, q, P; s, s' , s") + "2 cyclic"] 
- i [Gadik; T], s)Gbeiq; v\ s')Pcf{p] T]", s'Or^^y (k, q, p; s, s\ s" + "2 cyclic"] 
+ Gad{,k; T], s)Gbe{q] v', s')Gcf{p] r]", s")T^^^{k, q, p; s, s', s") \ , (56) 
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where "cyclic" refers to summation on the quantities obtained by permutations of 
(k, q, p) and, correspondingly, of (17, r/',//") and (s, s', s"). 

Since we are interested in the limit in which one of the external momenta vanishes, 
we notice that, due to the property of the trilinear vertex coming from translational 
invariance, limfc_!.o 7af)c(k, q, p) = 0, the only non vanishing contributions come from 
terms in which the vanishing momentum is associated to a </9-end of the trilinear vertices. 
This implies that the third line does not contribute, and that only terms in which the 
vanishing momentum is carried by a PS, and not by a propagator, have to be taken into 
account. 



Multiplying both sides by k ■ w (5_D(k), and using the Ward identities (47) and (51 ), 
we get 

j d^kk ■ w5d (k) Babe (k, g, |q + k\;ri, ri'ri") = j ds'ds"uaP^^^ {k) p • W X 
Gbe (g, r/', s') P^f {p, v", s") a:} (p, s', s") (e^" - e^') 

- Pbe (g, V, s') G,f {p, v", s") G-fl (g, s\ s') (e^' - e^") 

- Gbe (g, r/', s') G,/ (P, V\ s") (p, s'\ s') (e^" - e^') , (57) 

where we have used the fact that the PS carrying momentum k is in the linear regime. 
Only the first term in the first line and the first term in the second line of the expression 
within square brackets contribute to the final result: 

j d^k k-w 5Dik)Babe{k, g, |q+k|; r/, r,' , r]") = -q-w (e"' - e"") P'' {k)UaPbe{q; v' , v") , (58) 

The remaining terms in (57) cancel against each other, as can be seen from the relations 

j dsGae{q]ri,s)G~,^{q]s,7]') = 5ab5Div-v'), 
Pabiq; v, v') = Gae{q\ V, 0)G'M(g; 77', 0)ucUd P°(g) 

+ j dsds' Gae{q;V,s)Gbd{q;v\s')^cd{q;s,s') . (59) 
The identity (58) can also be expressed as 

\imBabeik,q,\q + k\;r],r]\7]") = -P\k)uaPbeiq;v' ,v")^ fe^' - e"") + 0(A;°) , (60) 

where by 0{k^) we indicate terms which do not have the q/k enhancement and which, 
in general, do not vanish for Vj' = t]" and are not proportional to the nonlinear PS 

Pbe{q;v',v")- 

In terms of the density contrast 6, and going back to conformal time, this relation 
rewrites 



lim Bs{k,q, |q + k|;r, t',t" 



-P,"(fc;r,r)P,(g;r',r") 



, ,,q-kD+(r')-/^+(r") 



A;2 



0{k^) 



(61) 
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where we recall that (r) is the growth factor at the time r. This relation is obtained 



by setting a = 6 = c= lin(60) and using eq. (14). Setting some of these indices to 2 
results in relations between correlators that include also the velocity field. 

We stress that this result is exact and nonperturbative, and merely descends from 
the requirement that the dynamics of the DM fluid is galilean invariant. In particular, it 
holds even in presence of vorticity and multi streaming, that is, when velocity dispersions 
and higher order moments cannot be neglected. Indeed, the derivation of the WI and of 



eq. (60) only requires that the action S and the initial PS in (29) are galilean invariant. 
As long as these conditions are satisfied, one can also introduce new degrees of freedom, 
such as velocity dispersion and vorticity, along the lines discussed in [22] , and integrate 
them over in the path integral without providing a source term for them. Therefore 
the range of scales to which this relation can be applied is in principle very large, being 
limited at very small scales only by the emergence of baryonic physics. 

5. Initial nongaussianity 



The relation (60) assumes that the initial fields are gaussian. It is interesting to study 
how it is modified in presence of an initial nongaussianity. For definiteness, we specify 
the current discussion to a nonvanishing initial bispectrum. Specifically, adding in the 



generating functional (29) the new term 
Z[Ja,K,;Vo] = 

j VipaVxh { ^ / ^'krf^q^aV |q + k|) Xa (k, 0) Xb (q, 0) Xc (-k - q, 0) + . . . | , 

(62) 

where dots denotes all the terms present in (29), we introduce a non vanishing x^ 



interaction in the free action, 

which, once attached to external propagators, communicates the effect of the primordial 
nongaussianity to the full PS, bispectrum, and so on. In particular, at zero-th order in 
labci we get the linearly evolved bispectrum of primordial origin, 

i^a^r (fc, |q + k|; r/, V, V') = gaM9he{v')9cf{ri") si (fc, |q + k|) . (64) 

In the following, we assume that, in the squeezed limit, the initial bispectrum 
factorizes a term proportional to the initial power spectrum of the non vanishing 
momentum, 

hm g, |q + k|) = fa (k) Pi (g; 0, 0) . (65) 

The initial nongaussianity is typically given in terms of a gauge invariant variable in the 
super horizon regime. For example, in terms of Bardeen's gauge invariant primordial 
gravitational potential we have 

lim (A;, g, |q + k|) = 4/nlP* {k) P$ (g) , (66) 
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where /nl = /nl*^' if the initial nongaussianity is precisely of the local form, while a 
different numerical factor is present for different other initial shapes. This relation is 
verified if the initial conditions are obtained from models of single field inflation [26] . 

In the linear regime, Bardeen's potential is related to our field (fa by the following 
relation, 

2k^T {k) 



(67) 



where the suffix at the denominator refers to the current time (where Hq = Hq, given 
that the scale factor is normalized to 1 today), and where T {k) is the transfer function 
of matter fluctuations, normalized such that T (fc) — )■ 1 when A; — > 0. Therefore, the 



relation (|64|) can be rewritten as a relation between the tree-level bispectrum and the 



linear PS's, 



limSr(^,g,|q + k| 



fc-s>0 



7],7] ,r] 



which is of the form (|65|) with 

fa{k) 



k^T (k) 



k'^T (k) 



(6J 



(69) 



(notice that Ua = 1). 



We are interested in how the relation (60) between the fully nonlinear bispectrum 



and PS's is modified at first order in the initial nongaussianity (given that the initial 
nongaussianity, if present, is small, higher order corrections will be very subdominant). 
Therefore, at the full nonlinear level in jabc, we should consider all the diagrams that 
contribute to the bispectrum that have at most one new vertex of the type given in 



eq. (63). The additional term in (62), involving equal-time fields, is Galilean invariant. 



and so it does not modify the Ward identities (47) and (51). Also the formal relations 



(56) and (59) are unmodified. The additional vertex however modifies the quantities 
Pab, Gab, ^ab, ^abc that enter in these relations. 




Figure 2. Diagrams contributing to the nonlinear bispectrum in the squeezed Umit, 
in which the top leg carries vanishing momentum. Only terms up to linear order 
in the initial nongaussianity are included. The external double lines denote the full 
renormalized PS and propagator, while the field in the upper vertex is still in the linear 
regime. 
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We discuss separately the two possible cases: either the three x ends of the 



additional interaction (63) are all linked to internal vertices in the loops, or one of 



them is directly connected to the external (k) through a propagator (where k is the 
smallest momentum in the bispectrum). In the first case, the contributions surviving 
the /c — )■ limit correspond to the first two diagrams in Fig. [2| which we denote by . 
The second case corresponds to the last two diagrams in the figure, which we denote by 
B^^ . The full bispectrum in the squeezed limit is then given by B^ + B^^ up to second 
and higher order corrections in the initial nongaussianity. 

Let us first discuss the B^ contribution. The first and second diagrams in Fig. u\ 



correspond, respectively, to those terms in the second and third lines of eq. (56) in which 



the vanishing k momentum is carried by an external PS, not by a propagator. As in the 



previous section, B^ does not receive contributions from the fourth line of eq. (56) and 



therefore it will be related to the nonlinear PS by the same identity (60) 

qk 



fe— ^0 



A;2 



(70) 



The difference between this expression and (60) is that now both the bispectrum and 



the full PS include also linear contributions in the initial nongaussianity. 

Let us now discuss the B^^ part. A generic diagram that contributes to it has the 
amplitude 

lim / i^a (k, T]) T^^l/ {k', q', | - q' - k'|) Xd (k', 0) Xe (q, 0) Xf (-k' - q , 0) . . . <^b (q, v') <^c (-k - 



lim>7a.(^) P {k)u. 



q,r/ 

P" {q', 0, 0) Xe (q', 0) Xf (-q', 0)---fb (q, v') V>c (q, v') 



In the first line, only the external fields and the B^x^ vertex have been written explicitly; 
to obtain the second line we have contracted (k) with one of the x from this 
interaction with a linear propagator, we have used ( [65| , and we have disregarded k 
when summed with q. The second expression is the amplitude of a diagram contributing 
to the full power spectrum Pbc (g; rj'rj") of the short modes at O (/nl) times a function 
(limfc_^o fi'aa' iv) fa' (^)) that depends only on the time and momentum of the soft mode 



(71) 




+ 




Figure 3. Diagrams contibuting to tiie full PS. 



An alternative way to see this is to realize that, in the squeezed limit of eq. (65) 



the new term of eq. (62) can be thought as a shift of the linear PS in the gaussian case 



of eq. (29) to 



P^M^PbMii-^fa{k)Xa{k,0)) , 



(72) 
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with k the vanishing mode. Therefore, all diagrams that contribute to B (at O^fMh)) 
can obtained by starting from a given diagram contributing to the full PS at O (/° 



and by replacing one of the linear PS that enter in this diagram by the second term at 



the RHS of eq. (72). The remaining Xa(k, 0) field is then connected to the external soft 
iPa field by a linear propagator. 

Therefore we get the nonlinear relation 



hm^ Bil {k, q, |q + k^, r/', r/") = f,{k) P,, (g; r/', r]") . 



fc-s>0 



(73) 



We stress that this is a relation between the full contribution to B^^ in the squeezed 
limit and the full power spectrum of the short modes. 



Once written in terms of the density contrast 6, the two relations (70) and (73) give 
the full bispectrum 



lim Bs {k, q, |q + k|; r, r', r") 

fc— s>0 

q ■ k (r') - jr") ^ e/NL^^.o^o' D+ (tq) 



fc2 



k^T{k) D+{t 

in terms of the full PS of the short modes, (where we have used the explicit expression 



Ps{k-r,T)Ps{q-,T\T") + 0{k\ /^l) , 

(74) 



for fa{k), eq. (69), and the fact that Ua = 1). 



6. Galieian invariance and perturbation theory 

Any computational scheme should reproduce the transformation properties of the 



correlators under a GT, eq. (39), in particular giving frame independent results for 



equal-time correlators. Since a GT can be seen as a coherent motion on very large 
(actually, infinite) scales, a non GT-invariant scheme would generally introduce a 
spurious dependence of small scale modes at momentum q on very large scale modes 
of momentum k <^ q, typically enhanced by powers of q/k over the true physical 
dependence. The role of the Ward identities is actually to guarantee that no such 
spurious effect are present in a given approximation scheme. 

In Sect. |4] we have seen how standard PT (up to 1-loop order) fulfills the Ward 
identities. In this section we want to investigate how the GT for correlators are 
obtained in standard PT and in particular, how the invariance is achieved for equal-time 
correlators. Then, in Sect. |8] we will discuss the same issues for some of the resummation 
methods proposed so far. 

As we have seen in Sect. |3} in a frame with non-zero average velocity, Vq, the field 
iPa has a non-vanishing expectation value - a tadpole in field theoretical language - given 



by ifa in eq. (34). This adds a new Feynman rule to the set of rules for the c.o.m. frame 
(Vo = 0) which, when properly taken into account, reproduces the transformation law 



of eq. (39). 



A convenient way to see this, is to rotate away the phases from the sources at 



the RHS of eq. (36), by performing the change of integration variables ifa{^^v) ~^ 



Galilean invariance and the consistency relation for the nonlinear squeezed bispectrum of large scale structw 
ipa{Kvy~'^'^°'^ , Xa(k,?7) Xa{KTl)e~'^-^°^ . It gives 

Z [J„ Kk, Vo] = e^/'^^'^^kj,(-M)^,(k,r,.Vo)gJ^>?rf='k ,^^(, [J„ K,- 0] . (75) 

The effect of the second exponential at the RHS can be seen by expanding in 
perturbation theory in the interaction vertex "^abc (the first exponential does not 
contribute to correlators with non vanishing external momenta), by using 

Z [Ja, Kf, 0] = e^^''^''^'^'^'^'P^"-('^''i'P)TOWWTOZo [Ja, K,- 0] , (76) 

where the free functional is obtained by setting '~^abc = in the path integral and 
performing the gaussian integral over ipa and Xa, [5] 

'1 



Zo[Ja, Kh] 0] = exp| — j d^lidridrj' 



+ iJa (k, 7]) gab (k; r], 7]') Kb (-k, r]') 



(77) 



The diagrammatic meaning of eqs. ( 75 ) and ( 76 ) is straightforward: starting from a given 



diagram for a given correlator in the Vq = frame and then inserting in all possible 
ways the tadpole given in Fig. |4] (weighting the insertion of n tadpoles by l/n\), the 
corresponding diagram for the Vq 7^ frame is obtained. 



k ; a 



-Vo • k II 6ab 



Figure 4. Inserting the tadpole in all possible ways in a diagram, integrating over 
the corresponding times, and summing over all diagrams contributing to a correlator, 
gives the GT of that correlator. 

Let US show how it works explicitly by considering the GT for the two-point 



correlator {i.e. the PS between unequal times), which, as we have shown in eq. (39), 
transforms as 

Pab (g; V, i) ^ e~^^-^^^''iv)-nv')] . (78) 

Let us first discuss the transformation of the linear PS, given in the Vq = frame by 

Pab (q; V, v') = 9ac{v)9MW)Pcdiq; o, o) = P'igWub , (79) 



where we have used eq. (32 ) and the property of the linear propagator, gabiv)^b = UaO{rf). 



Inserting the tadpole in the two possible ways, namely on the left and on the right legs, 
as in Fig. |5]we get 

- i q ■ Vo j dsdsT{s) [gac {v - s) (g; s, r]') - (g, r], s) gbc W - s)] 
= -2 q ■ Vo (T(r7) - T{r]')) P\q)UaUb , (80) 
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L 



+ 




Figure 5. One tadpole insertions on a linear power spectrum reproduce the first 



order term in the expansion in Vq of the linear part of ( 78 ) 



which reproduces the contribution to the transformation in eq. (78) at the hnear order 
in Vq. 

The O(Vo^) term is given by the three diagrams in Fig. [oj which give respectively 
the following three lines 

H q ■ Vo)' j dsdsT{s)ds'ds'T{s')S^gac {v " 9cd {s - s') (q; s' , r]') 

- 9ac iv - s) Pcd (q; s, s') gbd iv' - s) 
+ Pac s) gdc {s' - s) gbd iv' - s' 



Using the composition relation between the linear propagators 
9ab iv - s) 9bc (s - T]') = gac (v " v') 9{ri- s)e{s- r]') 



(81) 
(82) 



eq. ( |8l| rewrites 



;-^q-Vo) 



/•« rv rv 

ds I ds' — ds ds' 
Jo Jo Jo 



ds' / ds 
Jo 



d,T{s)d,>T{s') P\q)UaUb 



-ici-Yo{T{r^) - T{r,'))f P\q)UaUb. 



(83) 



The higher orders in Vq are reproduced by iterating the same procedure and using the 
property of the time integrals 



dsids-,T{si) / ds2ds2T{s2) 
Jo 



dS'ijids T (s^ J 



___L..t^ . . J- J 

Figure 6. Two tadpole insertions on a linear power spectrum reproduce the second 



order term in the expansion in Vq of the linear part of ( 78 ) 



The derivation above can be generalized beyond tree-level in perturbation theory. 
The crucial point is to realize how the tadpole insertions behave for closed loops. Indeed, 



by the same mechanism already seen in eq. (50 ), inserting in all possible ways the tadpole 



in the internal legs of a closed loop, and integrating over the time of the insertion, one 
obtains a contribution that is independent on the loop momentum, given by 



■ Vo T{si 



^5) 



Galilean invariance and the consistency relation for the nonlinear squeezed bispectrum of large scale structur 



where Pi and Si are the incoming momentum and the time argument of the external legs 
attached to the loop. Using this fact, one can realize that the contributions from the 
tadpoles factorize from a given loop diagram, provided they are attached in all possible 



ways to all the lines, and the transformation in eq. (78) (as that for any other correlator) 
is reproduced. Notice that this transformation holds for any diagram independently, that 
is, it is not necessary to sum all the diagrams contributing to the correlator at a given 



perturbative order to get the transformation in eq. (78) right. However, it does not 
mean that it is consistent, from the GI point of view, to discard some diagrams with 
respect to others of the same perturbative order, as we will discuss now. 

Consider again the two-tadpole insertions of Fig. |6} The three diagrams can be 
thought as the result of the following operation: take the 1-loop contribution to the PS, 
given in Fig. [7| and "open" one linear PS in all possible ways, namely, replace one linear 
PS, {ipa{q; s)ipb{p; s')) = (5z)(q + p)P°ft(g;s,s') with the product (^a(q; s)<^fe(p; s'), where 



<^a is given in eq. (34) This is indeed the effect of a change of frame, since, as we 
have seen, it corresponds to giving the expectation value to ipa- Put in this way, we 
realize that the O(Vo^) contribution of the GT of the tree-level PS can be obtained by 
taking the complete set of 1-loop diagrams and opening one linear PS. 

This can be generalized at all loops and at all orders in Vq. If we start from all 




Figure 7. One loop contributions to the power spectrum. 

possible diagrams with L loops that contribute to a given correlator, and we open P 
linear power spectra of that diagram in all possible ways, we end up with all possible 
diagrams with L — P loops contributing to that correlator, multiplied by the 0(Vo^^) 
term of the Taylor expansion of the GT phase, up to an overall combinatorial factor 
(2P)!!. This fact provides us a useful way to check the galilean invariance of a given 
approximation scheme. 

If we consider a GT-invariant correlator, such as the equal-time PS, then the sum 
of the opened diagrams should cancel. Physically, the cancellation reflects the fact 
that equal time correlators are insensitive to a global boost of particles velocities (see 



eq. (39) and next section). Therefore, locking the expansion in Vq to the PT expansion 



in 7af,c, as suggested by eqs. (75) and (76), provides the insightful link between GI and 
IR sensitivity, first elucidated in [20]. This link results in non-trivial relations between 
different contributions to a given quantity at a fixed PT order. In the 1-loop PS discussed 
before, the three diagrams of Fig. [7] exhibit a nontrivial cancellation between the mode- 

II Notice that, in the intermediate diagram of Fig. [rj it is possible to open two internal power spectra, 
obtaining the same result. This compensate the fact that this diagram enters with a relative 1/2 
multiplicity with respect to the other two diagrams in the Figure, while all the diagrams of Figure [6] 
have the same multiplicity. 
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mode coupling diagram in the middle and the other two diagrams, in which an external 
propagator is renormalized. 

If the opening procedure in a given set of diagrams for a galilean invariant correlator 
does not result in a complete cancellation, then that set is not insensitive to the effects 
of a GT transformation. 



7. Galielian invariance and IR sensitivity 



Lack of invariance under GT (in the locked Vo-7a6c expansion discussed at the end of 
the previous section) in a given calculational scheme can result in a spurious dependence 
of the short scale dynamics on large scales. To see this, we introduce a window 
function W {Xk) ~ 6' (1 — kX) (namely a regular function that however introduces a 
sharp difference between short and large wavelength modes), and define 

(k, 7^)^W (Xk) (k, r/) , (k, r^) = (k, r^) - (k, r^) , (86) 

and identically for the field Xa and for the sources Ja and Kh. For definiteness, we 
assume that we are in the frame of vanishing average velocity, Vq = 0. 

We are interested in the dynamics of the short wavelength modes in the background 



of the long wavelength modes. Under the separation (86), the cubic term —e^'jabcXa^b'^c 



then gives rise to interactions between these two types of modes. The interactions 
involving two long and one short wavelength modes vanish, due to the momentum 
conservation in the vertex. We obtain 



int,LS 



drie^'d^k.d-^^d-^qriabc (-k, -p, -q) 



Xa (k,^) </^b (P,r7) y^f (q,r7) 



+ Xa (k, ^b (P, ri) <^f (q, Tl) + 2xf (k, //) (p, v) (^f (q, r/) 



^7) 



The first term involves only long wavelength and approximately linear modes, and can 
therefore be neglected in this discussion, the second one is suppressed as A; p,q 
(since the vertex vanishes for A; — )■ 0). As for the third term, we assume that the long 
wavelength mode is still in the linear regime, and therefore we can express it as 



^b (p^v) 



-Hf 



■Ub 



-ie ^dfjTp ■ V (p) Ub ■ 



where V (p) is the Fourier transform of the physical velocity (corresponding to a peculiar 

proportionf 

2-5ac + 0(l) 



velocity TV(p)). The last term in (87) is then proportional to 

p ■ q 



7aic (-k, -p, -q)+7a2c (-k, -p, -q) 



2p2 



^^(k + p + q) , p<g.(89) 



The first term in this expression dominates once inserted into (87), giving 

^Sint,LS ^ - f dvd.T d'pd'q PqP;^(P) ^5 (_q _ p, ^) (q, r^) + . . . 



drid^c{ 



-d.Tq- I d^pYip) 



(90) 
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where we have written exphcitly only the contribution from the first term at the RHS 



of (89). However, once inserted in a correlator for the short wavelength modes, this 



expression has exactly the same effect as the tadpole term at the exponent of eq. (75), 
as it becomes even more clear if we make the replacement V(p) — )• Vo5_d(p)- Not 
surprisingly, from the point of view of short modes, the contribution of a velocity 
perturbation on ultra long modes is indistinguishable from that of a GT, and therefore 
it should have no effect on correlators that are galilean invariant. If GI is not properly 
taken into account, this spurious dependence of short scales on long scales can easily 
dominate over the physical one, induced by a long wavelength density perturbation as 
discussed, for instance, in [23], which, being proportional to 5d(p), and not to p5i?(p) 



as the velocity one, gives 0(1) terms from eq. (89). This would be the case, for instance, 
if, in the computation of the PS, one would neglect one of the diagrams in Fig. [7] and 
keep the other two. 

8. A galielian invariant resummation scheme 

As we have discussed in Sects. |6] and [7} standard PT respects GI: the computation 
of an equal time correlator at a finite PT order gives frame independent results once 
all the diagrams for that correlator at that order are taken into account. All different 
resummation methods, like RPT |3], the multipoint propagator expansion j^, closure 
equations [24], or the time-evolution equations of [17.j, amount to a reorganization of the 
perturbative expansion: in these schemes, a finite order approximation selects certain 
classes of PT contributions while neglecting others of the same PT order. Performing the 
computation in the new scheme at all orders would of course coincide with computing the 
result at infinite order in standard PT, and would therefore respect galilean invariance, 
but this is not automatically guaranteed when a finite order in the approximation is 
considered, as is mandatory on practical terms. 

Take for instance the computation of the PS in the RPT scheme discussed in [6]. 



The structure of the PS is the same as the exact expression in the second of eqs. (59), 
namely, 

Pabiq; V, v') = Gac{q; v, 0)Gbd{q; v', 0)ucUd P°(g) 

+ j dsds' Gac{q;r]iS)Gbd{q;v'^s')^cdiq;s,s') , (91) 

in which the full propagator Gab{k; ri,s) is replaced by an expression interpolating 
between 1-loop PT at small k and the large-fc result 

2 



G'^bik,V,v') =gabiv-v')ew 



(92) 



where = ^ f d^q ^ . The large-fc limit of the propagator emerges as the result of a 
summation, at any order in standard PT, of the so called 'chain diagrams' , i.e. diagrams 
obtained by inserting, in all possible ways, linear PS's carrying 'soft' momentum q <^k 



on the hard propagator line [1]. So, the first term in the expression at the RHS of (91) 
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is completely determined (modulo some arbitrariness in the interpolation procedure), 
while the second one, the mode-mode coupling, is given by loop diagrams in which the 
linear propagator is replaced everywhere by the non-linear one discussed above. Due to 
practical limitations, these corrections include at most some of the 2-loop diagrams in 
this new expansion scheme. 

Now, if we perform our 'PS-opening test', introduced at the end of Sect. |6]on the 
first term, we find contributions of the same sign, that do not cancel even in the equal 
time case, t] = rj . Opening, for instance, a single PS (in all possible ways) for all chain 
diagrams contributing to the resummed propagator gives an extra term 

- (k ■ Vo)^^ (T(r^) - T(0))^ {k, rj, 0) , (93) 
in the \aige-k limit, and therefore an extra contribution 

- (k ■ Vo)^ (T(r/) - T(0))^ GZ' {k, r^, 0) Glf {k, r^, 0) u,u,P'{k) (94) 



to the first term in (91) for the equal time PS. This result can be generalized 
to the opening of any number of linear PS, leading to the exponentiation of the 
— (k ■ Vo)^ ~ r(0))^ term. Now, positive contributions do actually come from 

opening PS's in the mode-mode coupling term at the second line of eq. (91)^ but they 
are not able to cancel the GT-induced term from the first line, such as eq. (94), as long 
as the mode-mode coupling therm is computed at a finite order in the new expansion 
scheme. Analogous considerations apply to the multi point propagator expansion [12]. 

The galilean invariance of nonperturbative methods based on the solution of 
evolution equations, such as [21] and [T7], requires an analysis of the RHS's of the 
differential equations. For instance, in the scheme proposed in [T7], the RHS of the 
time differential equation is GT invariant both in the — )■ and in the k ^ oo limit. 
However, the approximations done to treat the intermediate k range are not manifestly 
GT invariant. 

The numerical effect induced by mistreating galilean invariance should be carefully 
investigated case by case in each of these approaches, as each of them is based on 
a different reorganization of the perturbative expansion and on different truncation 
schemes. A naive estimate, obtained by computing the effect on the 1-loop PS of 
neglecting the contribution to the mode-mode coupling term coming from the IR modes 
responsible for the cancellation gives (at z = and for a ACDM PS) an effect around 
the percent at the upper extreme of the BAO range of k and growing as k"^ at higher fc's. 
Although this procedure likely overestimates the true effect, as the 1-loop IR cancellation 
is accounted for completely in all the methods mentioned above, we conclude that the 
galilean invariance of resummation schemes clearly deserves a careful consideration, 
especially if one aims to extend the range of validity of a given computational scheme 
beyond the BAO scales. 

% These contributions are those in which the flows of time from the two sides of the opened PS point 
towards the two different extrema at 77 and 77' of (91 ). PS of this kind are typically a few, compared to 
the infinite PS that can be obeyed in a resummed propagator. 
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In the remaining part of this section we discuss a resummation method introduced in 
[T7] . and show that is gahlean invariant. To our knowledge, it is the only computational 
scheme to pass the PS-opening test (and respecting Ward identities) among those 
discussed in the literature so far. A crucial role is played by counter terms, which 
should be properly taken into account at each order of the new expansion scheme to 
avoid over counting of the PT contributions. As we will see, they are also responsible 
for the fulfilling of the Ward identities derived in Sect. [4], and therefore, guarantee the 
galilean invariance of the expansion. 

The scheme can be considered as a variant of the RPT of |3l H]. The starting point 
is to define as the new tree-level propagator and power spectrum the quantities 



Gab (k; T], T]') = gab iv, V') exp 



(e^-e^Y 



9ab{V,V') 



ab ik;v,v'), 

dsds' gaM s)i::,'ik; s, s')Gtik; s', i) , (95) 



and 



GZ'ik; V, mn)Gt^{k- r,\ r,,^)P%k)uaU, 
+ / dsds'G^:^{k-i^,s)G'i:,\k-i,s')^^,\k-s,s') 



(96) 

for all values of k. These quantities are obtained by starting from the corresponding 
linear ones in standard PT, and by adding to them all possible chain diagrams, in which 
"soft" power spectra of momenta qi (namely, with qi -C k) are added to the hard line 
of momentum k. The approximation in which only these contributions are taken into 
account at the fully non-linear level is also referred to as the "eikonal" approximation 
[25] . due to its analogy with the resummation of soft gluon contributions in QCD. 
Therefore, we will refer to these tree level quantities as the eikonal propagator and PS, 
respectively, and, following [17] , we will indicate the new scheme as "eRPT" . 

To take these contributions into account, one can [17] add and subtract the 
quadratic expression 

iS, [$:^ {k- r/, r/') , {k- r/, r/')] = 
1 



drjdrj'd k 



--Xa (-k, r^) [k- r/, r/') Xb (k, V) - iXa (-k, r,) T.^^ [k- r/, r/') (k, r/' 



(97) 



in the exponent of the generating functional (29). The added term is included in the 



new "free" action, in such a way that, after after integrating out the <^a and Xa fields. 



Galilean invariance and the consistency relation for the nonlinear squeezed bispectrum of large scale structur 



the new free-path integral gives. 
Zf^'^lJa, K,-Yo] = expj - j d'kdr]dr]' 



+ tJa{Kr])G:^ik;v,v')K,{-Kv') 



+ I j d^kdr]J2 (-k, f]) ^2 (k, V, Vo) I . (98) 

The functions (k) and Jl'^^ (k) can be represented by diagrams in which the 
hard hne, carrying momentum of order k, is corrected by attaching to it soft power 
spectra in all possible ways so that the final diagram is IPI [T7]. The subtracted 



term, that is eq. (97) with flipped sign, is instead included in the new interaction term. 



The interaction lagrangian therefore consists of the cubic term from (23), plus the two 
counterterms —^at (k) (connecting two x) ^"^^ ^^ab^ (k) (connecting one and one x)- 



These counter terms avoid overcounting of the contributions already included in (98) 
By construction, eRPT and standard PT are equivalent at the infinite loop order. 



Repeating the same steps leading to (75) with the new splitting between free and 



interacting action, the generating functional can be expressed as 

\ /d,,r,'d3ke-*-^o(r(r,)-T{r,')) 



X 



<j,cik/j. I\ -Id -if> I ■■pcik/i l\ -Id -Id 

iK„{k,„)i./6(q.'))"c(p,r,) ^eik ^J^^ Qj _ 



X 



(99) 



Setting Vo = we get the new expansion in the eRPT scheme in the c.o.m frame. In a 



generic frame, comparing with (75), we see that there are new Vo 7^ 0-dependent terms. 



coming from the transformation of the counter terms. Therefore, a GT is realized by 
inserting, besides the tadpole term discussed in Sect. |6j the new contributions from the 
counter terms. 

This expression indicates how a correlator transforms under a boost of velocity Vo. 
We note that the transformation is realized by a series of mass insertions (identical to 
those in standard PT) and by phases on the counterterms. 



As an example, let us evaluate the GT of the propagator (95) up to first order in 
Vo. We obtain 

G:^ {k; r/, r/') ^ G^^ {k; r/, r/') - ^k ■ Vo /' dsdsT{s)G:^ {k; r^, s) Gf^ {k; s, r/') 



-zk-Vo / ds ds'G:^{k;7^,s)Kf{k;s,s')Gt{k-s',r^') [T (s) - T {s')] . 

Jri' Jri' 

(100) 

+ Integrating ipa out, now gives Xa (k, = drj'Jb (k, 77') Gl^^ (fc; ij' , 77); however, one does not need 
to change the frame fixing term of eq. (33), since S"'^**^ {k) vanishes at fc = 0. 
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Notice that for the eikonal propagators the composition property for hnear propagators, 
eq. (82), does not hold, and therefore we cannot take the product of propagators out of 
the integral on s at the first line, as we did in the PT case (see eq. (81 )). Here is where 
the role of counter terms in the second line is crucial. Indeed, we can write the second 
term at the RHS of the first line as 

- zk ■ Vo / ds 9,T(s)Gf (fc; 77, s) Gf {k- s, i) 



zk-Vo/ ds I ds' Gf^{k-r^,s)g~J{s,s')Gt{k-,s\i)[T{s)-T{s')] 



eik 



(101) 



where g^dis, s') = 6d{s — s') {Scdds' + ^cd)- Summing this expression with the second 
line of (|T00|) gives 

= tk-\o r ds f ds' (k; r/, s) G'^~\k; s, s')G'i^ (k; s', r/') [T (s) - T (s')] 

Jri' Jr]' 

= tk-Vo [ ' ds I ds' [Gf {k- 7], s) 6,bSD{s - v')Tis) 

-6adSD{v-s')T{s')G^^{k;s',v')], 
= -zk ■ Vo [T (rj) - T (V)] {k- V, v') , (102) 
which is exactly the linear term in the GT transformation for a propagator, according 
to the general rule of eq. (39). 

Proceeding in close analogy with the discussion for standard PT in Sect. |6} we 
can show that a cancellation between the Vo-dependent contributions takes place at 
any order in the eRPT expansion of the equal time PS, provided the counter terms are 
properly taken into account. 

We now discuss Ward identities in the eRPT scheme. The starting formal relations, 
eq. (45), as well as the derived ones as, for instance, eq. (46), are unchanged, since the 
generating functional from which they are derived is the same for PT and for eRPT. 
What is modified, on the other hand, is the way in which the GI enforced by Ward 
identities is implemented order by order. As in the calculation discussed above, the role 
of counter terms is crucial also in this context. 



For definitiveness, we consider again the Ward identity of eq. (46). At the tree 
level, the identity takes exactly the same form as for PT. Indeed the first term now gets 
two contributions, one from the inverse eikonal propagator and one from the counter 
term, whose sum gives back the inverse linear propagator g~^^, 

SD{p + q)G:r\q-v",v') 



r2peRPT 
" tree 



,/\ \'Pa=Xb=0 



+ 6d (p + q) Kti<i;v",v') 

tree i 



^^Pt 



\'fia=Xb=0 5 



Sxc (p, v") ^Vb (q, v'] 

while, since the tree-level effective action contains just one trilinear term, we get 



(103) 



SpcRPT 



6'r 



tree 



^Xc (p, V") (k, ri') 5(pb (q, T]') 6xc (p, v") (k, 5(pb (q, 



7Y- (104) 
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The one- loop level is less trivial. The next-to-lowest contributions to the second line of 
eq. (46) are given by diagrams such as those at the RHS of Fig. [ijwith the linear PS 



and propagator replaced by their full -eikonal- counterparts, while the vertices are kept 



at the tree level. As it was discussed after eq. (100), the two eikonal propagators (and, 
similarly, the eikonal propagator and the eikonal PS) meeting at the uppermost vertices 
do not combine into a single one, so that the LHS in Fig. [l] (in which the linear PS 
and propagator are also replaced by their full -eikonal- counterparts) is not reproduced. 



This is achieved, by the same mechanism of eq. (102), by adding to the LHS the three 
diagrams obtained by inserting a counter term in all possible ways, namely, at the 
left and at the right of the PS box and at the lower propagator. This mechanism works 
at higher orders, including higher loops and also the counter term. 

By similar arguments one can show that the eRPT scheme fulfills all the Ward 
identities imposed by galilean invariance, and is therefore safe from the point of view of 
spurious dependence on long wavelength velocity perturbations. 

9. Conclusions 

Galilean invariance is equivalent to the statement that long wavelength velocity 
perturbations completely decouple from short wavelength modes in the infinite 
wavelength limit. This is to be contrasted with the effect of long wavelength density (or 
velocity divergence) perturbations, which can be seen as a modification (renormalization) 
of the background evolution for the short modes, as discussed in [23]. Therefore, 
computing a given correlator in a frame different from the c.o.m. one, in which 
the average velocity of particles is Vq 7^ (in practice, by considering the tadpole 
contributions of Sect. [6]) should give the same result as the one obtained by 'gauging 



away' Vq by a GT (and then considering the phase factors of eq. 39). This line of 
thinking leads to the Ward identities obtained in this paper, and is the counterpart in 
the LSS context of what is usually done for gauge symmetries in quantum field theory of 
for conformal invariance in inflation, where the used of Ward identities and consistency 
relations is well established EZl EHl [211 EOl ED E2] • 



Nonlinear relations such as eq. (74) can be of use in many respects. From the 
computational point of view, as discussed in Sect. [8| they provide consistency checks 
for nonperturbative approaches aiming at extending the range of validity of PT towards 
smaller scales. In this paper we have identified a well defined expansion scheme, eRPT, 
which, thanks to the role of counter terms fully passes the test. The numerical feasibility 
and performances of eRPT, as well as the implementation of explicitly GT invariant 
approximation schemes to the approach in [T7j, will be analyzed in a future work [33] . 

We stress again that, since these relations descend merely from the GI of the DM 
fluid, they hold even at scales in which the single stream approximation is broken, and 
therefore apply also to more refined schemes, in which velocity dispersion is taken into 
account, such as [221, and to the construction of effective field theories, as in I34|. 



Last, but not least, the relation between the nonlinear PS and the nonlinear 
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bispectrum in the squeezed limit, expressed in eq. (74) should be eventually probed 
in numerical simulations, and in real observations, once large redshift surveys will 
be completed and the correlators at different times will be measured. The different 
time behavior of the two dominating terms in (74) can possibly allow to extract the 
contribution from the primordial nongaussianity, as opposed to that induced by the 
nonlinear dynamics during the formation of the structures. A quantitative study of this, 
and of other relations that follow from GI (for instance, between higher point correlators, 
or between correlators including the velocity field) is definitely an interesting topic for 
future research. 
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